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1. Introduction 

Rapid advances in experimental techniques of laser cooling and trapping made it possible 
to achieve Bose-Einstein condensation (BEC) in weakly interacting systems. One of the 
many interesting aspects of BEC in dilute gases is the dynamics of growth and decay 
of the condensate PlElEj- Bose-Einstein condensates can have long lifetimes ranging 
from 2 s 4J up to more than 10 s E]. This finite lifetime is mainly caused by 
inelastic collisions between condensate atoms and by collisions with particles from the 
background gas, resulting in atom loss and heating of the system. Many experiments 
with condensates can be performed on a time scale short compared to the condensate 
lifetime, such that decay effects are not important. This is, however, not always the case 
and some experiments specifically focus on condensate decay [IJ. Hence, it is important 
to model the decay of a condensate in detail, taking into account all processes which 
significantly contribute. In literature condensate decay due to inelastic collisions and 
collisions with background gas atoms as well as heating effects are discussed. To our 
knowledge, the transfer of atoms between the condensate and the thermal cloud occuring 
dynamically as a consequence of thermalization of the system has not yet been taken 
into account. We will show that this process plays an important role when the number of 
atoms in the condensate is of the same order or smaller than the number in the thermal 
cloud. Its effect could only be neglected in experiments with large condensate fractions 
which are often obtained by removal of most of the thermal part with an rf-knifc, or in 
situations in which our assumption of fast thermalization is not fulfilled. 

The present work concentrates on the decay of condensates in the presence of 
a considerably large thermal fraction. We assume thermal equilibrium during the 
decay of the condensate which is justified in most experiments with BEC in dilute 
atomic gases. Since elastic collision rates are typically large (> 10^ s~^ IH IHl), 
one expects thermalization to occur very rapidly jZj compared to the rate of change 
in thermodynamic variables during the decay. Using the simple condensate growth 
equation j2] we investigated numerically the effect of dynamical disturbances caused 
by atom loss and heating from inelastic collisions on a system originally in thermal 
equilibrium. We found that the system stays very close to thermal equilibrium under 
conditions typically encountered in the experiments with dilute atomic gases. 

Although our model incorporates two- and three-body collisions it does not take 
into account subsequent secondary effects. For instance, the effects of induced local 
variations in the mean-field interparticle interaction jH] are neglected. Also secondary 
collisions of reaction products with atoms of the condensate or the thermal cloud are 
not accounted for. Avalanches, recently discussed by Schuster et al are therefore 
beyond the scope of this work. They may be incorporated later. 

This letter is organized as follows. In Sec. El we derive a simple analytical formula 
for the decay of a condensate of non-interacting bosons including only losses due to 
collisions with background gas particles. The expressions derived in this section show the 
existence of transfer of atoms between condensate and thermal cloud and the significant 
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reduction in condensate lifetime that may occur. In Sec. |21we present a more complete 
dynamical model of condensate decay using the mean-field theory for weakly-interacting 
bosons and including also losses by inelastic two- and three-body collisions. In Sec. IH 
we present, as an example, numerical simulations of condensate decay in metastable 
helium. 



2. Atomic transfer 



Apart from elastic collisions that keep the system in thermal equilibrium, in this section 
we assume the atoms to undergo only collisions with background gas particles, inducing 
atom losses. For simplicity we consider a system of non-interacting bosons. At a 
temperature T below the critical temperature Tc = hio[N/g3{l)Y^^ /k the number of 
atoms Nj- in the thermal cloud is given by PU] 

iv.=rf)(g'. (1) 

where oo = {uixUJyUJz is the averaged frequency of the harmonic potential trapping the 
atoms, denotes the number of atoms and gn{u) = XlfcLi '"'^/^") dsi^) ~ 1-202. The 
energy Et contained in the thermal cloud is JU] 

with a = 7r'^hujg^'^^^{l)/30. To describe the equilibrium state we choose instead of the 
temperature T and the total number of atoms N as dynamical variables the number of 
atoms in the condensate Nc and the number of atoms in the thermal cloud A^^- 

For the background gas particles we assume room-temperature energies and a 
uniform distribution over the volume occupied by the trapped atoms. As the average 
kinetic energy of a background gas particle is much larger than the energy of an atom 
in the trap, the collision cross section does not depend on the latter and each collision 
removes one atom from the trap. After such a collision the total energy of the system 
E is depleted on average by the mean energy per atom E/N. With these assumptions 
we obtain for the loss rate 

N = --N = Nc + Nt = --{Nc + Nt), (3) 

T T 

with T denoting the trap lifetime. For the rate of change in the total energy E of the 
system accompanying the atom loss we have 

1 • • 1 

E = — E = Eq + Et = — {Ec + Et), (4) 
r r 

where we used the energy contained in the condensate Eq = EqNc with Eq = 
^h{iJx + ujy + ujz) being the ground state energy of the trap. Using Eq. ^ the time 
derivatives Et = ^aN^^NT and Ec = ^o^c can be substituted in Eq. (jH) and in 
combination with Eq. © we then find 

AT 1 /vr 1 - s^Nt/Et 3 

Nt = — Nt — ; —, — ^ Nt, (5) 

T 4/3 - EqNt/Et 4r ^' ^ ^ 
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where we used EqNt/Et 1, neglecting the energy of a condensate atom compared to 
the average energy per atom in the thermal cloud. Substituting in Eq. Q yields 

Nc = -^(^Nc + ^NTy (6) 

This simple analysis shows that for a condensate coexisting in thermal equilibrium with 
a thermal cloud in a trap neither Nc = —Nc/t nor Nt = —Nt/t holds for the decay 
induced by collisions with background particles. Conservation of energy and number 
of atoms combined with rapid thermalization inevitably results in transfer of atoms 
from condensate to thermal cloud thereby enhancing the decay rate of the condensate. 
Especially with a considerable fraction of thermal atoms in the system this affects the 
decay of the condensate. Only in the limit of a large condensate fraction Eq. Q becomes 
Nc = -Nc/t. 



3. Interacting model 



In order to parametrize the equilibrium state of the gas, we will use the temperature T 
and the number of atoms in the condensate Nc as independent variables fully describing 
the state of the system. We will start with the stationary "two-gas" model proposed 
by Dodd et al ^l], in which atoms of the thermal cloud do not affect the condensate 
described by the stationary Gross-Pitaevskii (GP) equation. The thermal cloud atoms 
do not interact with each other (except for thermalization) but they are influenced by 
the condensate through the mean-field potential 2t/onc(r). Here, the contact potential 
Uq = Airh^a/m expresses the binary interaction between atoms with scattering length 
a and mass m, and nc(r) denotes the spatial density of the condensate which can 
be calculated from the GP equation. We additionally simplify the two-gas model by 
describing the thermal cloud semiclassically, replacing discrete states when evaluating 
statistical averages by a continuum of states. This way we obtain analytical formulas 
for the atomic density ^^(r) and energy density exir) in the thermal cloud: 

d^p 1 



{27rhy 



93/2 



exp 



2m 



+ K//(r)-/i AT 



,-{Veffir)-t,)/kT 



(7) 



d^p 



exp 



-kT\ 5(5/2 I e 



V,ffir)-f,)/kT , ^ 



(8) 

where A = ^y2'n-h?/mkT is the thermal de Broglie wavelength, /i is the chemical 
potential and Ve//(r) = V{r) + 2Uonc{r). The explicit form of nc{r) and /i follows 
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from the GP equation. They depend on the number of atoms in the condensate Nq- 
By integrating the densities nc and nx over spatial degrees of freedom we obtain the 
number of atoms in the thermal cloud Nt{Nc,T) and the energy of the thermal cloud 
Et{Nc,T). 

We are interested in the time dependence of Nc and T. As before, we will first 
consider the dynamics of the total number of trapped atoms and total energy E a.s a. 
function of Nq and T. The advantage of this approach is that we do not need to state 
transfer terms explicitly; the transfer will follow from our analysis automatically. 

The total loss rates and E are related to A^^ and T via 



dT \dN, 



■ OEt- fdEr 



(9) 

with /i = dEc/dNc- The reparametrization is straightforward as we have already found 
explicit expressions of almost all terms appearing in Eq. 0. 

In order to describe the dynamics of the system we must consider all relevant 
processes affecting its state. We include three main dynamical effects that may 
cause losses of atoms from the system: two-body inelastic collisions, three-body 
recombinations and collisions with background gas. In the following we will neglect 
all secondary collisions. Then the total atomic loss rate A^ is given by [T^ : 

- N = + 2x j d^r |^^n^ + 2ncnT + + 

3^ / ^''^ (^^^ + l^c'^r + S^c'^T + 4) , (10) 

where x ^ are two- and three-body collision rate constants and r is the lifetime of 
the trap. 

Each term in Eq. ()10|) corresponds to a loss process that may occur in the system. 
For example the term ^^n^(r)nT(r) expresses the probability density of a three-body 
recombination between two condensate and one thermal cloud atom to take place at 
point r. In this event we lose all three atoms from the system, hence the factor 3 in 
front of the integral in Eq. The energy lost in this particular process consists of two 
terms. First, an energy 2/i carried by the two lost condensate atoms. Second, the energy 
carried by the thermal cloud atom just before the recombination. We assume that the 
change in internal energy per lost thermal particle is equal to the average energy per 
atom eT{r)/nT{r) at point r. Then the energy loss rate by this process becomes 



(11) 

Analogously, one can derive the proper rates corresponding to the other possible 
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collisions. Finally we end up with the expression for the total energy loss rate: 



E = \et + fiNc) +2x j dV 



eT t n \ 

— [Uj, + ncrtT} + 



+ 



3^ / d\ 



— I ^^c'^T + 2ncnrp + 1 + 

^ 1 + j^J^c^T + ricriT^ . (12) 
In equations (fTUj) . and (fT^ A^, i?, A^t, Et, and are functions of Nc and T only. 

4. Results and conclusions 

In this section we present some numerical solutions of the set of equations ©, (UHl), and 
(|12|). We have performed simulations of the condensate decay for the conditions of the 
Paris experiment on metastable helium an initial number of atoms in the condensate 
of A''c(0) = 5 X 10^, a trap lifetime r = 35 s, s-wave scattering length a = 16 nm, trap 
parameters uj^ = = 1090 Hz and uo^ = 115 Hz, two-body inelastic collision rate 
X = 1.5 X 10~^^ cm^/s |13j, and three-body recombination rate ^ = 4 x 10"^"^ cm^/s 
|14j . The initial number in the thermal cloud is either Nt{0) = 5 x 10^ (Fig- ^ or 
Nt{0) = 2 X 10® (Fig. 12)), with initial temperatures of 1.5 /iK and 2.5 /xK, respectively. 
This gives densities still outside the regime where avalanches occur. The numerical 
solution Nc{t) is shown together with results for a much simpler approach in which the 
transfer is absent. For the loss rates Nc and Nt again Eq. (|TO|l is used, setting to zero 
the thermal or condensate density, respectively. Comparing both curves with the same 
initial conditions it follows that, when the existence of transfer is neglected, the lifetime 
of the trap is overestimated by 50% when the initial fraction of atoms in the thermal 
cloud is 0.5 and by more than a factor of 2 when this fraction is 0.8. The third (dotted) 
line in each figure shows the decay in the absence of a thermal cloud. The lifetime 
increases by a factor of 2.5 in Fig. [T]and a factor of 8 in Fig. |21 due to the absence of 
inelastic collisions between condensate and thermal cloud atoms. In both models the 
decay rate is larger with a thermal cloud, but in the full model the effect is enhanced 
by atomic transfer. 

Another feature of condensate decay that can be studied on the basis of our model 
is the time dependence of the temperature during the decay. Usually it is assumed that 
the temperature is constant p. In both figures an inset shows the numerical solution 
for T{t). Indeed, the temperature remains constant within 4%. 

We have presented an equilibrium model of the decay of a Bose-Einstein condensate. 
Our analysis has shown that the assumption of sustained thermal equilibrium leads to 
transfer of atoms from the condensate to the thermal cloud which can significantly 




Figure 1. Decay of a condensate of metastable helium atoms with. Nc{0) = 5 x 10^ 
applying our model (solid line) in comparison with a simpler model neglecting 
thermalization (dashed line). The initial number of thermal atoms Nt{Q) = 5 x 10^ 
[r(0) w 1.5 /iK]. The dotted line represents the decay of a pure condensate: iV7-(0) = 
(T = K). Inset: dynamics of the temperature in our model. 
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Figure 2. Decay of a condensate with a large thermal cloud: A'j'(O) — 2 x 10^ 
[T(0) « 2.5 ^K]. Further details as in the caption of Fig. ^ 



enhance the condensate decay rate. This effect could be seen by an experimental 
examination of the decay rate as a function of the fraction of thermal atoms. 
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